Prom closed to open ID Anderson model: Transport versus spectral statistics 
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We show that the transport properties of a one-dimensional Anderson model of finite size with 
weak disorder can be effectively expressed in terms of the repulsion parameter < /3 < oo in the 
level spacing distribution of eigenvalues of the corresponding closed system. This result stems from 
the detailed numerical analysis demonstrating that the normalized localization length of eigenstates 
is nothing but the parameter /3. We give the analytical expressions for the mean transmission 
coefficient (T) and its variance Var(T), as well as for (InT) for any value of /3 and degree n of 
coupling to continuum. The numerical data fully correspond to the analytical predictions. 
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Introduction. In spite of a remarkable progress in un- 
derstanding statistical properties of quantum systems, ei- 
ther deterministic or disordered, one of the important 
problems still waits for the detailed analysis. The spe- 
cific question is: to what extent one can predict scat- 
tering properties of a complex open system, if we know 
global properties of eigenstates and energy spectra of the 
corresponding closed system? This problem was solved 
for a specific case of closed systems with maximal chaotic 
properties described by fully random matrices. A proper 
mathematical tool in these studies is based on effective 
non-Hermitian Hamiltonians of a certain structure [1|, 
and many analytical results have been obtained, see, for 
example, [2|-y]. The key point in this method is that 
the scattering matrix of an open system is expressed in 
terms of eigenvalues and eigenfunctions of the related 
closed system along with their decay amplitudes. 

A much more difficult problem emerges if the closed 
system is not fully chaotic being characterized by an addi- 
tional parameter related to the degree of chaos. Recently, 
this problem was analyzed numerically in Refs. [5| where 
the global characteristics of scattering or signal transmis- 
sion through a system have been studied in dependence 
on two control parameters, the degree of internal chaos 
and the strength of coupling to continuum. In particular, 
it was found that, independently of the degree of chaos, 
the increasing continuum coupling leads the system from 
the quasistable regime of isolated narrow resonances to 
the "super-radiant" regime of overlapping resonances co- 
existing with long-lived compound states. However, the 
specific features of this evolution may critically depend 
on the degree of chaos and therefore the observation of 
the signal transmission provides important information 
on regular or chaotic character of intrinsic dynamics. 

The above studies have been performed with the use of 
random matrices (canonical Gaussian ensembles or two- 
body random interactions). Some statistical assumptions 
may be questionable in application to realistic physical 



systems. Below we study the ID Anderson model, pay- 
ing main attention to the relation between the scattering 
properties of an open model and those of eigenstates and 
spectral statistics of the closed model. We discover un- 
expected effects that give a new insight to the problem 
of scattering through finite disordered systems. 

The model. The tight-binding Anderson model is often 
used to describe electron propagation in random media. 
In the ID case the corresponding Hamiltonian with di- 
agonal disorder takes the standard form. 



Hn 



iy{5. 



rn,n -t-1 



+ <5, 



771, n — 1 



(1) 



Here v is the hopping amplitude connecting the nearest 
sites (in what follows we set i^ = 1); the site energies e^ 
are assumed to be uniformly distributed in the interval 
—W/2 to W/2 giving rise to the disorder variance a^ = 
W^ /12. Our interest is in the transmission properties 
through the samples of finite size N with the arbitrary 
coupling amplitudes \p^ and VT^ connecting the left 
and right edges with attached semi-infinite ideal leads (in 
which e„ = 0). In the case of zero disorder open tight- 
binding models have been studied in Refs. [a, |7|. As for 
non-zero disorder, so far, the main interest was related 
to the statistics and distributions of resonances for one 
open channel [8|, |9[ ; the relation of the resonances to the 
transport properties was studied in [10|. 

Without disorder, the spectrum of the closed chain 
consists of Bloch waves with the nodes at the ends and 
energies inside the band \E\ < 2. In the limit iV — >■ oo 
and for weak disorder, cr^ <C 1, all eigenstates are expo- 
nentially localized with the characteristic length loo{E) 
given by the Thouless relation |ll| . 
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This expression is valid everywhere apart from the vicin- 
ity of the band edges, \E\ =2, and band center, E ^ 0. 
Level repulsion in a closed model. To quantify the de- 
gree of chaos in the finite samples with no continuum con- 



plings, 7^ = 7^ = 0, we employ the well known results 
of random matrix theory (RMT) . The maximal chaos in 
deterministic quantum models with chaotic behavior in 
its classical counterpart is characterized by the Wigner- 
Dyson (WD) distribution P{s) of normalized spacings s 
between the nearest energy levels. In the opposite case 
of integrable classical counterparts the level spacing dis- 
tribution is typically close to the Poisson distribution. 
Therefore, one can take the distribution P{s) as a mea- 
sure of chaos in the closed Anderson model. 

For zero disorder, the energy spectrum near the band 
center, e„ = 0, is equidistant, so that P{s) -^ 6{s — 1), 
the eigenstates are extended and regular. In the limit 
of strong disorder, all eigenstates are effectively localized 
on the scale of the sample size iV ^ 1, thus resulting 
in the Poisson distribution. In between these limits, for 
l^ « A'^, one can expect that the eigenstates are both 
extended and "chaotic" . By the last term we mean that 
the components of eigenstates are uncorrelated, with the 
distribution close to the Gaussian. In this region P{s) 
is expected to be close to the WD-distribution. Thus, 
when passing from zero to strong disorder, the repulsion 
parameter /3 in P{s) changes from 13 = oo to f3 = with 
/3 w 1 in the intermediate case. 

In order to describe the evolution of P{s) in a large 
range of the values of /3, we use, instead of the Brody 
interpolation [l^ , the expression suggested in Ref. [l3| , 
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where /(/3) = ^'^21^ (l ^ f) " 0.16874. Here z = Trs/2 
and the parameters Bi and B2 are determined by the nor- 
malization conditions, L Pp{s)ds = L sPi3{s)ds = 1. 
The function /(/?) is constructed in such a way that for 
the values /3 — 1;2; 4, corresponding to the Gaussian en- 
sembles of random matrices of a given symmetry (orthog- 
onal, unitary and symplectic, respectively) it is close to 
the expressions for P{s) obtained in the RMT for those 
ensembles [12[. As shown in Ref. 14 1, for these values 
of /3, the dependence ^ is more accurate (in a whole 
range of s) than the WD surmise typically used in the 
literature. Note that for /? = Eq. ^ coincides with 
the Poisson distribution, and for /3 = cx) it reproduces 
the delta-function. With Eq. ([3]) we performed a numer- 
ical study by changing the degree of disorder for a closed 
chain in a large range of the control parameter x ~ loc/N 
with loo defined by Eq. ([2]), see examples in Fig. [1] 

Our conclusion is that Eq. (J3|) gives an amazingly good 
correspondence (supported by the x^ criteria) with the 
numerical data in a very large range of x. It should be 
noted that Eq. ([3]) was obtained in Refs. [13|,llJ] by using 
the analytical expressions by Dyson |15| derived for the 
classical model of two-dimensional charged particles on 
a ring with the temperature 1//3. For specific values of 
the inverse temperature /3 = 1,2,4 the partition function 
giving the probability to find the particles at specific po- 




FIG. 1: Examples of P{s) for 0.02 < \E\ < 0.2 (excluding 
the energies very close to 0), with A'' = 1000. The data are 
obtained for 120 disorder realizations with the x^-fit to Eq. Q 
for r = 40 bins. Confidence levels are given in parenthesis. 



sitions on the ring in this Coulomb gas model is the same 
as for the canonical Gaussian ensembles. Therefore, one 
can make an unexpected conclusion that the distribution 
P{s) for the Anderson model on a finite scale A^ can be 
described by the Dyson Coulomb gas model for which s 
is the distance between the nearest particles on a ring. 

Localization length vs. repulsion parameter. Now we 
can establish the relation between x and /3. It is qual- 
itatively clear that they express the same phenomenon 
of gradual transformation of standing waves into local- 
ized states. All the data are summarized in Fig. [2] We 
see a precise linear dependence between x and /3 in a 
whole range of x values independently of the chosen en- 
ergy range. The fit of the data as /3 = Ax + C gives the 
slope A = 2.3 ± 0.1 with C essentially zero. This result, 
obtained carefully with the use of x^ statistical criteria, 
shows that the repulsion parameter (3 is nothing but the 
properly normalized localization length /qo. 
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The factor 2.34 in Eq. (j4]) can be attributed to the fiuc- 
tuations of components of eigenstates. 
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FIG. 2: Repulsion parameter /3 versus x = loc/N for _B ~ 0, 
(circles), and _E ~ — 1 (squares), see details in Fig. [T] 

The localization length loo can be defined due to the 
Shannon entropy, S = ~^n=i^^ri^'^Wn with w„ = -0^^. 
For completely random states with the Gaussian distri- 
bution of i'n one gets, S = ln(A^/2.07). Thus, it is 



convenient to introduce the normalized entropic local- 
ization length, (In = 2.07 N exp{S), where (...) repre- 
sent an ensemble average. With this definition we have 
djv = N for fully chaotic eigenfunctions occupying N 
sites. Our numerical data show that the onset of strong 
chaos occurs when cIn « 2.1loo UM which is equivalent 
to /? « cLn /N K, 1. Therefore, one arrives at the same re- 
sult through an analysis of the spectrum as through the 
eigenstates. 

Open model: Non-Hermitian Hamiltonian. The scat- 
tering properties of open systems can be formulated ll | 
with the effective non-Hermitian Hamiltonian [y, |8|, |lO|, 
3, where V{E) = E/2 - {i/2)y/4~~W, 
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This exact expression is valid for any disorder e„, contin- 
uum coupling 7^,7^, and energy E. Near the center of 
the band (the general case of any energy, — 2 < £■ < 2, is 
studied in Ref. [16[) it can be written as 
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where W{E) is defined by the coupling amplitudes, 

AL,R^E) ^ (^^■«/.)V2[i _ i,V4]^/^(4a^ + 4^). (7) 

Here Hmn from Eq. ([1]) describes the internal dynamics, 
while the non-Hermitian part W{E) is factorized in terms 
of the coupling amplitudes A'^{E) between the internal 
states \n) and open decay channels labeled by c = L, R. 
The non-Hermitian Hamiltonian ([6]) allows us to con- 
struct the scattering matrix S in the space of channels. 
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where the reaction matrix K is defined as 
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and ipm{Ej) is the m-component of the eigenstate \j) of 
the closed Hermitian Hamiltonian ([1]). 

In the case of weak disorder, cr^ ^ 1, the strength 
of coupling to the leads is characterized by the coupling 
parameter, k^ = 2tt^'^/ND, calculated via the average 
scattering matrix. For this definition the channel trans- 
mission coefficient r"^ = 1 — | (5''^'^) p is maximal for perfect 
coupling when k'^ — 1. Here D is the mean level spac- 
ing at the center of the energy band in a closed chain, 
k"^ = 0. Below we consider the equiprobable coupling, 
7*^ = 7, k"^ = K. The average scattering matrix can be 
written in the form (5^"^) = (1 — k)/(1 -I- k) and vanishes 
for perfect coupling. In all numerical simulations we take 
A'^ = 1000 sites and combine an ensemble average over a 



large number of realizations with a spectral average over 
1000 energies across an energy window \E\ < 0.2. 

Transport characteristics. It is known (see, for exam- 
ple, [l7|) that, for continuous weak random potentials in 
ID geometry, all transport characteristics for finite sam- 
ples of length TV depend on the ratio N/Loo [E) . Here Loo 
is the back-scattering length of scattering states that can 
be associated with the localization length ^oo of eigen- 
states in infinite samples. This is the core of the single 
parameter scaling (SPS) [l8|. Unlike the models with 
continuous potentials, the SPS is questionable for the 
tight-binding model ([T|). It is shown that at least at the 
band center, E = Q, the SPS fails ^). 

Below we demonstrate that the main results obtained 
for ID disordered models can be also applied to the An- 
derson model provided the energy E is not very close to 
band center. For example, in the case of perfect cou- 
pling, the average transmission coefficient and its second 
moment (g = 1 and 2, respectively) are given by [17 1 
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which together define the variance of the transmission, 
Var(T) EE (r2) - (T)2. Here x = loo/N with loo defined 
by Eq. ©, and (...) stands for the ensemble average. A 
simpler relation emerges for the self-averaged logarithm, 
(InT) = ~2N/lao. Quite often, this expression serves as 
the definition of the localization length /qq. Our numer- 
ical data for (T), (T^), Var(r) and (InT) manifest excel- 
lent correspondence with Eqs. ([TO)) . 

For non-perfect coupling we derived the following for- 
mula for (InT) which is valid for any coupling k, 
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Our data nicely correspond to this expression for 0.05 < 
K < 10 and 0.2 < /? < 22 with /3 defined by Eq. (4); 
/3 = 1 marks the strongest chaos in a closed system. 

Another quantity of theoretical and experimental in- 
terest is the reflectivity factor F , 
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where R and T are the reflection and transmission co- 
efficients, respectively, and the term (5'^^)^ takes into 
account direct scattering emerging for a non-perfect cou- 
pling [5[. For perfect coupling, k = 1, we have (5''^'^) = 0, 
and Eq. ([T^ with R = 1 — T can be written explicitly 
with the use of Eq. PH)) . For non-perfect coupling, k ^ 1, 
we found an approximate relation. 
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that gives quite a good description of the data for oq = 
2.8, fli = 1.4, see Fig. [3l For /3 — >■ oo the expression 
F = (1 — k)^/(1 + kY can be confirmed analyticahy. Our 
data demonstrate that the reflectivity factor strongly de- 
pends on the degree of internal disorder. It is remarkable 
that regardless of the coupling there is a sharp change in 
the value i^ at /3 « 1 corresponding to the onset of strong 
"chaos" in the closed system. Thus, measuring experi- 
mentally F one can observe the transition from extended 
to localized states on the scale of the sample size N . 
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FIG. 3: Reflectivity factor versus tlie parameter /3 = 
2.34Zoo/A^ for A'^ = 1000 and different couplings. Each sym- 
bol is obtained by an average over 1000 disorder realizations. 
The inset shows a closer view of the indicated section with a 
logarithmic scale on the vertical axis. 

Conclusion. We studied the transport properties of the 
ID Anderson model in dependence on the degree of in- 
ternal chaos and strength of coupling to continuum. We 
found that the level spacing distribution P{s) for a closed 
model of finite size N is well described by the phenomeno- 
logical expression ([3]) in which the repulsion parameter 
P changes from /? = to /3 = oo. This expression is 
originated from the two-dimensional Coulomb gas model 
with the temperature 1//3, and gives the distribution of 
spacings between nearest charged particles moving on a 
ring. This fact may be used for further analytical studies 
of the spectrum statistics in the finite Anderson model. 

In the closed model we established the linear relation 
between the repulsion parameter /3 and normalized local- 
ization length, loo/N, of the eigenfunctions. Thus, the 
repulsion parameter (3 reflects the transformation of the 
wave functions from localized states to extended stand- 
ing waves. In a narrow region with /? '^ 1 the properties 
of the model can be adequately described by the RMT. 

Opening the system at the ends we used the effective 
non-Hermitian Hamiltonian to study the transport 
properties of the model. We demonstrated that the 
average transmission coefficient and its variance entirely 
depend on the parameter /3 for a given strength of 
continuum coupling. Finally, the reflectivity factor 
(the ratio of average fluctuative reflection to average 
transmission) is shown to be a universal function of 
the same parameter. We conclude that the trans- 
mission properties of the Anderson model provide a 



reliable tool for getting information on the degree of 
randomness and spectrum statistics in the closed system. 
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